Spin-dependent Klein tunneling in graphene: Role of Rashba spin-orbit coupling 
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Within an effective Dirac tlieory the low-energy dispersions of monolayer graphene in the presence of Rashba 
spin-orbit coupling and spin-degenerate bilayer graphene are described by formally identical expressions. We 
explore implications of this correspondence for transport by choosing chiral tunneling through pn and pnp 
junctions as a concrete example. A real-space Green's function formalism based on a tight-binding model is 
adopted to perform the ballistic transport calculations, which cover and confirm previous theoretical results 
based on the Dirac theory. Chiral tunneling in monolayer graphene in the presence of Rashba coupling is 
shown to indeed behave like in bilayer graphene. Combined effects of a forbidden normal transmission and 
spin separation are observed within the single-band n "H> p transmission regime. The former comes from real- 
spin conservation, in analogy with pseudospin conservation in bilayer graphene, while the latter arises from the 
intrinsic spin-Hall mechanism of the Rashba coupling. 

PACS numbers: 72.80. Vp, 72.25.-b, 73.23.-b, 73.40.Gk 
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I. INTRODUCTION 

After the first successful isolation of monolayer graphene 
(MLG) was announced,' intriguing properties based on its 
low-energy excitation that mimics massless, gapless, and chi- 
ral Dirac fermions were intensively investigated.^'^ Spin-orbit 
coupling (SOC), on the other hand, is the key ingredient of 
semiconductor spintronics"*-^ that was undergoing a rapid de- 
velopment before the rise of graphene.^ The question about 
the role of SOC effects in graphene then naturally emerged, 
including the proposal of graphene as a topological insulator,^ 
which attracted the attention of various first-principles-based 
studies.**"'" 

SOC in MLG includes an intrinsic and an extrinsic term. 
The former reflects the inherent asymmetry of electron hop- 
ping between next nearest neighbors^ (i.e., a generalization 
of Haldane's model"). The latter is induced by the elec- 
tric field perpendicular to the graphene plane, which can be 
externally controlled, and resembles the Rashba model'^ '^ 
for the two-dimensional electron gas. Agreement has been 
achieved, based on first-principles calculations,^'"' that the in- 
trinsic SOC term opens a gap of the order of 2/1/ sa 24 /zeV, 
while the Rashba SOC removes the spin degeneracy and cre- 
ates a spin-splitting 2Xr at the K and K' points that has a lin- 
ear dependence on an external electric field E with the slope 
of about lOOjLieV per V/ A of £. Under a strong gate voltage, 
the Rashba coupling may in principle dominate the intrinsic 
SOC in MLG.'* '" 

The low-energy spectrum of MLG plus the Rashba cou- 
pling (MLGh-R) was derived by Rashba,'^ based on the Kane- 
Mele model^ ( i.e., an effective Dirac Hamiltonian). An ear- 
lier work by one of us'^ started with a tight-binding model 
(TBM) and obtained an equivalent form of the low-energy 
expansion,'^ 



with |5k| fl ^ 1. Recall for comparison the low-energy spec- 
trum of bilayer graphene (BLG),^'" 



£mlg+r (q) « At 5 [\/ (3f«)^ + (3ffl -qf + v {3tR)] , (1) 

which also agrees with expressions given in Refs. 9 and 10 
when A/ = 0. Here /i, v = ±1 are band indices, t and Ir are 
nearest-neighbor kinetic and Rashba hopping parameters, re- 
spectively, fl w 1 .42 A is the bonding length, and q = K+8k 



£blg {q)^^l^^{^Jrli+{^ta■qf + VYl), 



(2) 



where Yi is the nearest-neighbor hopping between the two 
graphene layers. Note that the next nearest-neighbor inter- 
layer hoppings 73 and 74 do not influence the band dispersion 
near K. The completely different mechanisms of (i) pseu- 
dospin coupling between carriers from the two graphene lay- 
ers of BLG through interlayer hopping Ji and (ii) real-spin 
coupling between up and down spins within MLG through 
Rashba hopping Ir happen to lead to an identical mathemati- 
cal form in Eqs. (1) and (2) that can be clearly mapped onto 
each other'^''^ with Ji ^ ^tR as sketched in Fig. 1. This un- 
ambiguously implies that low-energy physics in MLGh-R and 
BLG should behave similarly. 

In this paper we tackle the question of whether the trans- 
port in MLGh-R behaves as in BLG by choosing the issue 
of Klein tunneling^'^'^"'^' (or, in general, chiral tunneling) as 
a concrete example. Chiral tunneling in graphene has been 
shown to exhibit completely different behavior in MLG and 
BLG based on the Dirac theory.^^ Tunneling at normal inci- 
dence in MLG shows a suppression of backscattering, which 
resembles the original Klein paradox in relativistic quantum 
electrodynamics^^ and hence the name Klein tunneling, while 
in BLG it shows a perfect reflection, which is strictly speaking 




FIG. 1. (Color online) Schematic of the pseudospin coupling through 
7i in BLG (left panel) and real-spin coupling through Ir in MLG 
(right), which lead to an identical low-energy dispersion near K and 
K'. 
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a consequence of forbidden interband transition also due to the 
chiral nature of graphene. The theoretical discussion of chiral 
tunneling so far focuses mainly on spin-independent tunnel- 
ing through pn and junctions, while SOC effects 
are less discussed. In addition, the relevant theoretical 
understanding so far is based on Dirac theory, which is valid 
only for the Fermi level close to the charge neutrality point and 
allows only to consider certain relatively simple systems. A 
recent study discussing the interplay between the Aharanov- 
Bohm effect and Klein tunneling in graphene, started with a 
TBM,^'* but the nanoribbon type of the leads used in that work 
may have edge effects included that can be very different from 
the bulk properties of graphene. A more transparent theoreti- 
cal study of chiral tunneling in graphene directly bridging the 
analytical Dirac theory and the numerical TBM computation 
is so far missing and deserves consideration. 

In the present work, we re-treat this issue of chiral tunneling 
in graphene based on the TBM and show a unified description, 
allowing for a broad range of geometries and complementing 
the existing results based on the Dirac theory. Straightfor- 
ward generalization to the case of MLGh-R reveals a spin- 
dependent tunneling behavior in close analogy with that in 
BLG, with the role of pseudospin in BLG replaced by real 
spin in MLGh-R. Specifically, a combined behavior of spin- 
Hall-based spin separation and suppression of normal trans- 
mission will be shown. 

This paper is organized as follows. In Sec. II we briefly 
summarize the theoretical formalism applied in the present 
calculation, namely, real-space Green's function formalism in 
noninteracting bulk graphene. In Sec. Ill we show our TBM 
results including the consistency with the Dirac theory, a di- 
rect comparison between BLG and MLGh-R, and a deeper dis- 
cussion of the MLGh-R case. We review also briefly the re- 
cent experimental progress on the Rashba spin splitting and 
Klein tunneling in graphene in Sec. IV, and finally conclude 
in Sec. V. 



II. FORMALISM 

A. Tight-binding model for "bulk" graphene 

We choose the TBM for describing the electronic prop- 
erties of graphene, which is a well established way to treat 
graphene numerically. For spin-degenerate MLG, the Hamil- 
tonian reads 



in Eq. (3) but will not be considered in the present work due 
to the minor role it plays in the bulk transport properties for 
low-energy excitation. 

Spin-orbit interactions can be incorporated into the TBM 
by altering the spin-dependent hopping between nearest and 
next-nearest neighbors,^'^^ modifying Eq. (3) as 

^LG+R = Y.V,a^c}c, + £ [-ta° + itR{a x d,;),] cj. 

{ij) 

(4) 

Here a" is the 2x2 identity matrix, Ir is the Rashba spin- 
orbit hopping parameter, d,; is the unit vector pointing from 
site j to /, and a = (c'^, (7^,(7^) is the vector of (real-) spin 
Pauli matrices. We take into account only the extrinsic SOC 
and neglect the intrinsic term in order to highlight the role of 
the Rashba SOC. 

For spin-degenerate BLG, we consider 

Mlg= L ^<:!5-7iE(^L«ij + H.c.), (5) 

"1=1,2 j 

where J^j^^q is J^/^mlg given by Eq. (3) of the mth graphene 
layer, a„,j {b,nj) annihilates an electron on sublattice A (B) in 
layer m ~ 1 , 2 at unit cell j (that contains two sublattice sites 
belonging to A and B), and the interlayer coupling strength 
7i w 0.4eV corresponds to the nearest neighbor hopping be- 
tween the two MLG layers. Further interlayer hopping terms, ^ 

-r4lj(«2j«i > + ^L-^iJ+ H-c.) and -r3l;(4j^i >+ H-c), 
are not considered in the present calculation, since they do not 
influence the low-energy excitation. Throughout the presen- 
tation of the numerical results in Sec. Ill, the kinetic hopping 
parameters will be fixed at f = 3eV and 7j = 0.39 eV, while 
the value of the Rashba hopping parameter Ir depends on the 
context. 

For the simulation of bulk graphene, we impose the Bloch 
theorem along the transverse direction with periodicity W. 
This is equivalent to considering a nanoribbon and modify- 
ing the hopping between atomic sites connected through the 
periodic boundary conditions by a Bloch phase factor e''^^^ 
with a Bloch momentum kg,^^ as schematically shown for 
MLG in Fig. 2. At the same time the Bloch momentum 
is the component of the electron's momentum perpendicu- 
lar to the nanoribbon, hence defining the propagation angle 
(j) ~ sin^^kB/kp), where kf is the Fermi wave vector To be 
consistent with the literature related to Klein tunneling based 
on the Dirac theory, in Sec. Ill we will refer to the Bloch mo- 
mentum as ky. 



=^LG=EV,c!c, -f £c!cy, (3) 

('■J> 

where the operator (c,) creates (annihilates) an electron at 
site / (including both sublattices A and B). The first sum in 
Eq. (3) runs over all the atomic sites in the considered region 
with on-site potential Vj, and the second sum runs over all the 
pairs of neighboring atomic orbitals with kinetic hop- 
ping parameter t (w 3eV). The next nearest neighbor kinetic 
hopping term, usually characterized by t' « O.lf, can be added 



pik^W ^ik^W ^ik,,\V ^ikyW ^ikyW ^tkyW ^ikyW ^ikyW ^JkyW ^ikyW ^ikyW 




L lead S region R lead 



FIG. 2. Schematic of a minimum tight-binding model that simulates 
a bulk MLG up to nearest neighbor hoppings with W = 3a. Fur- 
ther nearest neighbor hoppings can be accounted for by enlarging 
the transverse periodicity W to at least 6a. 
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In the present calculations, we will apply a minimal TBM 
by imposing the periodic boundary conditions on a zigzag 
nanoribbon with chain number A'^ ~ 2, that is, periodicity of 
W = 3a (as the case sketched in Fig. 2). The present model ap- 
plies equally well for metallic armchair ribbon (chain number 
Na being a multiple of 3) with periodic boundary conditions, 
but the minimal model would require Na = 3 (i.e., periodicity 
ofW = 3V3a). 



B. Brief summary of real-space Green's function formalism 

We consider open systems connected to the outer world by 
two leads (see Fig. 2). According to the real-space Green's 
function formalism^^ we numerically calculate the Green's 
functions of our system, 

Gf = [E-Hs-r^"±iri]-\ (6) 

where the self-energies of the leads (E'/" = L'^^" + T.'^l") re- 
flect the fact that our system is open. The powerful recipe 
constructed in Ref. 36 for graphene handles a lead as a semi- 
infinite repetition of unit cells and allows for incorporating 
any kind of lattice structure and one-body interaction such as 
SOCs. The transmission probability for an electron traveling 
from lead L to lead R is given by the Fisher-Lee relation^^-^^ 

7>L = Tr(riGjrKFs"), (7) 

where the trace is done with respect to the lattice sites. The 
spectral matrix functions Ti/^ are given by the lead self- 
energies as Fi/fi = /(E[/„ - E^/^). 

For a given Bloch momentum ky and a given Fermi energy 
Ep [subject to a Fermi wave vector kp via Eq. (1) for MLGh-R 
or Eq. (2) for BLG], the incoming propagation angle (p of the 
electron wave can be defined as ^ — sin^^ {ky/kf). The angle- 
dependent transmission function T {(j)) is obtained from Eq. 
(7), which can be generalized to a spin-resolved version. 

III. TRANSPORT RESULTS 

In this section we present numerical results of our tight- 
binding transport calculations. We first show the consistency 
of our tight-binding calculations with the existing effective 
Dirac theory in Sec. Ill A. A direct comparison between BLG 
and MLGh-R will then be shown in Sec. IIIB. Finally, Sec. 
IIIC is devoted to MLGh-R for pn junctions, in particular the 
role of Rashba SOC for chiral tunneling. 

A. Consistency witli Dirac tlieory 

We first consider tunneling in graphene without SOC and 
confirm existing results, limited to low energy excitations, by 
our tight-binding calculations. We pick two pioneering theo- 
retical works to demonstrate the consistency explicitly. Con- 
sistency with recent works of tunneling in graphene hetero- 
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FIG. 3. (Color online) Tunneling through a barrier for (a), (b) 
MLG with Ef = 3lakf /2 = 81.6meV and (c), (d) BLG with Ef = 
{3takf/2)-/Yi = IT.lmeV. In (b), red (light gray) and blue (dark 
gray) curves correspond to Vq = 196.8meV and Vq = 280.3meV, 
respectively. In (d), red (light gray) and blue (dark gray) curves cor- 
respond to Vq = 48.7meV and Vq = lOO.TmeV, respectively. In both 
cases the barrier width is D = lOOnm and the incoming Fermi wave 
vector is kf = 2;r/50nm^', as considered in Ref. 22. 

junctions in the presence of SOC'^"'' has also been checked, 
but is not explicitly shown here. 



1. Chiral tunneling in MLG vs BLG 

Tunneling in MLG and BLG behaves quite differently 
as mentioned in Sec. I and pointed out by Katsnelson et 
al?^ For a quantitative comparison we consider a barrier of 
width D = lOOnm and the incoming Fermi wave vector kf ~ 
27r/50nm-' as in Ref. 22 for both MLG and BLG [see Figs. 
3(a) and 3(c)]. Note that in order to exactly match the bar- 
rier width, we set the bonding length a — (4-\/3)^' nm, which 
differs from the realistic value of about 1. 42 A by only less 
than 2%, so that the number of hexagons used here amounts 
toD/(\/3fl) = 4x [D]nm = 400. 

The resulting transmission probabilities as a function of the 
incident angle are depicted in Figs. 3(b) and 3(d). They 
reproduce the results of Fig. 2 in Ref. 22 almost perfectly, if 
we choose slightly different Ef and Vq, to which the trans- 
missions at finite angles are sensitive. The remaining tiny 
difference between our TBM results and their Dirac theory 
results^^ simply reflects the basic difference between the two 
approaches: For graphene the effective Dirac theory is vahd 
only for energies close to the Dirac point, while the TBM is 
suitable for the entire energy range. 

Note that the maximal values of the transmission functions 
in Fig. 3 are 2, since the valley degeneracy is automatically 
incorporated in the tight-binding formalism. Later when we 
take spin also into account, the maximum of the transmission 
function will be 4. The transmission probabilities calculated 
by the Dirac theory always have their maximum of 1 due to 



4 




2 


— Eq. (8) 






---TBM 




1.5 






S 1 










\ 


0.5 















-90 -60 -30 30 60 90 

^ (deg) 

(a) 



-90 -60 -30 30 60 90 

4> (deg) 



(b) 



FIG. 4. (Color online) Klein tunneling in MLG through a pn junc- 
tion with a (a) sharp and (b) smooth interface, (a) Comparison be- 
tween TBM (dashed line) and Eq. (8) [solid green (gray)] showing 
perfect agreement (Ef = 80meV). (b) Comparison between TBM 
(long and short dashed) and Eq. (11) [solid green (light gray) and 
red (dark gray)] for kfd ^ 6.16 (Ef = 200meV) and kpd ^ 1.54 
[Ef = 50meV), respectively. 



the normalized incoming wave, unless a proper degeneracy 
factor is taken into account. 



2. Klein tunneling in MLG: Sharp vs smooth interface 

Tunneling in MLG through a pn junction exhibits probabil- 
ity one at normal incidence and is called Klein tunneling. In 
experiments, a graphene pn junction can be realized by using a 
backgate, which tunes the carrier density (and hence the Fermi 
level) globally, and a topgate that tunes locally the carrier den- 
sity, equivalent to the potential step Vq at the other side."*^ The 
carrier densities on the two sides can be controlled to be of 
opposite signs, forming the pyi junction. In between, however, 
the variation of the carrier density is never abrupt in reality. 
Cheianov and Fal'ko showed, based on the Dirac theory, that 
the interface of the /injunction actually matters.^"* They con- 
sidered symmetric pn junctions (i.e., Vq = 2Ef) with sharp 
and linearly smooth interfaces, which we briefly review and 
compare with our TBM results in the following. 

a. Sharp interface For a symmetric pn junction with a 
sharp interface [see the schematic in Fig. 4(a)], the transmis- 
sion probability as a function of was written as^^ 



r(0) = cos^0, 



(8) 



which does not depend on the potential step height. This sur- 
prisingly simple expression matches our TBM result always 
perfectly as long as Vb = 2Ef, as shown in Fig. 4(a). 

For a step potential with arbitrary height Vb ^ 2Ef, the 
transmission probability as a function of the incident angle 
^ and the outgoing angle Q can be derived as 



2 cos cos 9 

1 +cos(0 + 0) 



(9) 



which agrees with our TBM calculation equally well as the 
symmetric case (not shown). The two angles and Q are 



connected to each other due to conservation of transverse mo- 
mentum by 



sin d 



\Ef-Vq\ 



sin^, 



(10) 



where i = for nn' or pp' and —1 for np or pn. Equation 
(9) clearly recovers the symmetric pn junction case of Eq. (8) 
when choosing s = — I and Vo = 2Ef in Eq. (10). Note that 
in the case of \Ef — Vol < I^fI- the Fermi wave vector in the 
outgoing region is shorter than that in the incoming region, 
and an additional constraint for has to be applied to ensure 
|sin0| < 1 [i.e., (/) < with(^), = sin-'(|£f -yo|/|£F|)]. 

Previously it has been stated that the single-valley Dirac 
picture, based on which Eqs. (8) and (9) are derived, is not 
equivalent to the TBM."*' The difference in their work, how- 
ever, becomes noticeable only when the distance between one 
of the involved energies and the Dirac point exceeds roughly 
300 meV. In our simulation, indeed the deviation for the sym- 
metric pn junction case with, say Ef = 300 meV, is less than 
0.5%. The agreement of our TBM and the Dirac theory there- 
fore confirms that the intervalley scattering, which is mainly 
responsible for the nonequivalence at high energies, is indeed 
negligible. 

b. Smooth interface For symmetric pi^ junctions with a 
linearly varying region of width d [see the schematic in Fig. 
4(b)], the analytical derivation for the transmission probability 
within the Dirac theory yields^^ 



r(^) =exp 



-K Sin 

2 ^ 



(11) 



for kfd ^ 1.^^ This formula, together with the validity cri- 
terion kpd 1, are tested by our tight-binding calculations 
shown in Fig. 4(b), where two sets of parameters are consid- 
ered. For kfd K, 6.16 we find very good agreement with Eq. 
(11), while the result for kfd « 1.54 exhibits noticeable devi- 
ations from the analytical prediction at large angles |0|. The 
smoothing function was assumed in their work as linear but 
the reality might be much more complicated, which is then 
not accessible by the Dirac theory but again straightforward 
by our tight-binding calculation. Nevertheless, the exponen- 
tial form of Eq. (11) is still a good description regardless of 
the actual form of the smoothing function, as we have numer- 
ically checked. What really matters is only the product kfd. 

Unlike the sharp pi^ interface, a compact form of transmis- 
sion probability for the asymmetric case does not exist so far. 



B. pnp junction: BLG vs MLG+R 

We next show the direct correspondence between BLG and 
MLGh-R by considering exactly the same potential barrier 
and incident Fermi energy as in Fig. 3(d) for BLG, and set 
3f« = 7i = 0.39 eV for MLGh-R here. (A discussion with 
weaker, realistic tR will be continued in the next section.) The 
total transmission shown in Fig. 5 for MLGh-R indeed resem- 
bles the curves in Fig. 3(d) for BLG, as expected due to the 
identical form of their low-energy dispersions (1) and (2). The 
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most important feature of chiral tunneling in BLG, forbidden 
normal transmission, now appears also in the case of MLG+R. 
In BLG, T {(p =0) =0 was understood as the consequence of 
pseudospin conservation. For MLG+R, T {(p =0) =0 can be 
expected as the consequence of real-spin conservation. In- 
deed, this can be demonstrated by computing the nonequilib- 
rium local spin density, which can be obtained from the lesser 
Green's function,'*^ considering two cases, < Ef < 3tR and 
—3tR < Ep < Q, both with = 0. Within this single-band 
transmission, the local spin densities for positive and negative 
Ef point to opposite directions, indicating that normal inci- 
dence transmission between n and p regions will be forbid- 
den. 

Next we discuss the spin-resolved transmission. The quan- 
tization axis is chosen as the out-of-plane direction, so that 
the transmission of, for example, Ti^^- means the probability 
of an incoming +S- electron ending up as an outgoing —S- 
one. Since the incoming angle dependence (j) of the transmis- 
sion probabilities are analyzed, we define 7| = 



- as the 



t,; Too' = Tf + Ti = T is 



transmission ability of the electron (or t spin), and vice 
versa. (Alternatively, one can also analyze the outgoing angle 
dependence and define as + 7||, not used here. Either 
way, the total transmission Y,a.a 
ensured.) 

The choice of quantization axis z is not necessary but facil- 
itates relating the present spin-dependent tunneling in MLG 
with the issue of intrinsic spin-Hall effect previously dis- 
cussed in semiconductors. The spin-resolved transmission 
curves shown in Fig. 5 exhibit opposite lateral preference of 
the t and | electron spins, which is an intrinsic spin-Hall 
mechanism due to the Rashba SOC. In a semiconductor two- 
dimensional electron gas (i.e., a continuous system rather than 
discrete as in the TBM), such an intrinsic spin-Hall deflection 
of opposite electrons can be easily explained by the con- 
cept of a spin-orbit force based on the Heisenberg equation of 
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FIG. 6. (Color online) Transmission T at normal incidence (ky = 
0) as a function of potential step height Vo for tunneling through a 
pn junction in MLG+R. The leftmost solid band diagram above the 
main panel corresponds to the incoming n side. The five ticks on the 
Vo axis correspond to the above five dashed band diagrams for the 
outgoing side. 
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(12) 



Here Jif ~ p^/2m + {aR/h){pya'' — p.,<7-^ ) is the continuous 
two-dimensional Hamiltonian with Rashba SOC, r and p are 
the position and momentum operators, Ur is the Rashba cou- 
pling parameter (rather than the hopping one, f^), and C7' is the 
sign of the 5, spin component. The 7| and 7| curves shown in 
Fig. 5 therefore reveal a combined effect of forbidden normal 
transmission due to conservation of real spin and the intrinsic 
spin-Hall deflection that can be understood by Eq. (12). 

A few remarks are due before we move on. To connect 
BLG with MLG+R we put 3tR = Yi = 0.39 eV, which is ap- 
parently far from reality. In general the Rashba splitting in- 
duced by electrical gating is roughly of or less than the order 
of 100 jiieV (see Sec. IV). Fermi energy lying within this split- 
ting, which is also our main interest, projects to a much shorter 
Fermi wave vector kp, leading to a much longer up to a few 
or a few tens of microns in order for kpd 1 to be valid. This 
implies that the influence of the interface on the tunneling in 
MLG+R is normally negligible, unless d is that long. In addi- 
tion, tunneling through a pnp junction will also require a long 
barrier width D for electrons subject to such a short kp', other- 
wise, the barrier is merely a weak perturbation to the electron 
due to its long Fermi wave length. Based on these remarks, we 
will focus in the next section only on pn junctions in MLG+R 
with a reasonable Rashba hopping parameter. 



FIG. 5. (Color online) (a) Angle-resolved total transmission T for 
tunneling through a pnp junction in MLG+R with the same barrier 
height Vo, barrier width D, and Fermi energy Ef as used in Fig. 3(d) 
for BLG, and a substitution 3/r = = 0.39eV. (b) and (c) show 
spin-resolved transmission probabilities for Vq = 48.7meV and Vo = 
100.7 meV, respectively. 



C. pn junction in MLG+R 

In the following we demonstrate in detail the role of Rashba 
SOC in tunneling through a potential step in MLG+R. The 
Rashba hopping parameter will be fixed to tR ~ 30 /ieV and 
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the Fermi energy in most cases to Ep = 2tR, which lies within 
the spin-orbit spHtting (see Fig. 1). 



1. Normal incidence 

We begin with the case of normal incidence, ky = 0. In Sec. 
TUB we have discussed the one-band transmission selection 
rule (i.e., n ^ p transmission is forbidden). The transmission 
from the left side at Fermi energy < Ep < 3f^ to the right 
side with potential Vo is expected to be zero whenever a single- 
band « — > /? transmission is attempted. Indeed, as shown in 
Fig. 6, a zero transmission gap of T as a function of Vb is 
found. The gap lies in the interval of Ef <Vq < Ef + 3tR, cor- 
responding to the single-band n ^ p transmission. Note that 
contrary to the valley-valve effect in zigzag nanoribbons,'*^^^ 
the gap shown here arises solely due to a bulk property. 



2. Angle- and spin-resolved transmission 

We proceed with angle- and spin-resolved transmission and 
consider first the trivial case with Ef ~ 0.5 meV well above 
the Rashba splitting 3ffl = 90/ieV, as shown in Fig. 7. In 
this case the maximum of T = + is 4 since two spin 
subbands and two valleys are involved in transport. The to- 
tal transmission curve resembles the expected cos^ ^ behavior 
as discussed in Sec. Ill A 2, showing that the Rashba effect 
plays only a minor role. The spin-resolved 7^ and curves 
differ only slightly at |0| = sm-\k'^/kf') k. 56°, where k'^ 
and k°p^ are the inner and outer radius of the two concentric 
Fermi circles, respectively. Tunneling in BLG with Ef well 
above 7; behaves similarly (i.e., the interlayer coupling 7; in 
BLG no longer plays an important role in the process of chiral 
tunneling when the transport occurs at Ef ^ 7; ), as we have 
numerically checked. In other words, the chiral tunneling in 
BLG with Ef > and in MLGh-R with Ef > 3f^ recovers 
the Klein tunneling behavior as in MLG. 

Of particular interest is the nontrivial case with \Ef\ < 3?^. 
As a test, we first consider Vo = as shown in Fig. 8(a). In the 
absence of the potential step, the total transmission function 
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FIG. 7. (Color online) Angular dependence of total (T) and spin- 
resolved (T^ and T^) transmissions for Ef = 0.5 meV well above the 
Rashba splitting 3tR = 90/ieV. 
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FIG. 8. (Color online) Angular dependence of total and spin-resolved 
transmissions through a pn junction in MLG+R with (a) zero poten- 
tial and (b) finite potential. Parameters used are given above. 



T reaches its maximum of 2 (one spin subband times valley 
degeneracy of two) for any angle 0, as it should. The opposite 
lateral deflection tendency of the '[ and | spins is again clearly 
seen and can be explained based on Eq. (12) as discussed in 
Sec. IIIB. 

The most important case is that of Fermi energy Ef E 
(0,3?^) and potential height Vb G {Ef,Ef-\-3tR). A specific 
example with Vo ~ 100 /ieV is shown in Fig. 8(b), which ex- 
hibits the combined effect of the forbidden normal transmis- 
sion [r(0 = 0) = 0] and spin-Hall deflection. The number 
of high transmission peaks is always two.^^ Compared to the 
previous trivial case (Ef > 3tR, Fig. 7) where and Ti^ do 
not significantly differ, the separation of the opposite f and 
I spins is distinctly enhanced. Whether this could be a new 
type of intrinsic spin-Hall mechanism in graphene deserves a 
further investigation, and is left as a possible future direction. 

We summarize the discussion of angle- and spin-resolved 
transmission by mapping r(^,Vo) in Fig. 9. Four different 
transport regimes can be identified: 

1 . Vo < 0, single n band to single/multiple n band(s) trans- 
mission regime. 

2. < Vo < Ef, single n band to single n band trans- 
mission regime; distinct spin-resolved and T^, 
and high total T limited by a critical angle 0^. = 

sin-i(|£f-Vo|/|£f|). 

3. Ef < Vo < Ef + 3tR, single n band to single p band 
transmission regime; combined effects of forbidden 
normal transmission and spin-Hall deflection. 

4. V(,> Ef + 3tR, single n band to multiple p bands trans- 
mission regime. 
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FIG. 9. (Color online) Transmission through a pn junction in MLG+R as a function of incident angle (j) and potential step height Vq- Four 
transmission regimes can be distinguished: (i) Vq < 0, (ii) < Vq < Ep, (iii) Ef<Vo<Ef + St^, and (iv) Vq > Ep + 31^, with Ef = 60 fleY 
and 3fR = 90jUeV. 



Note that a vertical scan in Fig. 9 at ^ = corresponds to 
Fig. 6, and horizontal scans at Vq = and Vq = 100 /ieV to 
Figs. 8(a) and 8(b), respectively. These four regimes will be 
helpful in the following discussion of conductance. 



3. Integrated conductance 

Finally, we calculate the conductance of the pn junction in 
MLG+R by integrating T {(p), or equivalently, T (ky), with re- 
spect to the transverse Bloch momentum, 

111 r'^F 

G=-^/ T{ky)dky, (13) 

2kF J-kp 

where the prefactor ensures the maximal value of the 
Landauer-Biittiker type ballistic conductance to be /h times 
the maximal number of modes. We compare the conduc- 
tance of the pn junction in MLG [Ir = 0) and in MLG+R [tR = 
30 /xe V) as a function of the potential step height Vq, as shown 
in Fig. 10. Since the Fermi level is fixed to Ep = 60/xeV for 
both cases, the transport for Ir = will involve two spin and 



Ep Ef + 3tR 




Vo (MeV) 



FIG. 10. (Color online) Integrated conductance of the ballistic pn 
junction in MLG with = and MLG+R with tR = 30/xeV. 



two valley degeneracies, leading to the maximal G of 4-e^/h, 
while in the case of tR — 30 /ze V only one spin subband is pro- 
jected, leading to the maximal G of 2e^/h. The maximal G 
occurs always at Vq = that corresponds to an ungated clean 
bulk graphene. Zero conductance, on the other hand, occurs 
at Vq = Ep since no states at the outgoing region are available 
at this charge neutrality point. 

Different transmission regimes can be distinguished based 
on our previous discussion for Fig. 9. For Vq G [0,60]/ieV 
(« — !• n transmission), the rise of Vq shrinks the Fermi circle at 
the outgoing region and hence introduces a critical transverse 
momentum, outside which the transmission is suppressed due 
to the lack of out-going states. The critical transverse mo- 
mentum reduces linearly with Vq for MLG due to the linear 
dispersion. The conductance G, Eq. (13), therefore reduces 
also linearly with Vq. In the presence of the Rashba SOC, the 
low-energy dispersion becomes quadratic, and so does the re- 
duction of G with Vo in MLG+R. 

For Vo G [60, 150] /leV (« — > p transmission), the conduc- 
tance of MLG rises faster than that of MLG+R, possibly due 
to the help of Klein tunneling. At Vq = 150/ieV, a sudden 
jump (or a shoulder) occurs in the case of MLG+R since the 
second spin subband at the outgoing region starts to partici- 
pate in transport. This jump does not occur in the MLG case 
since both spin subbands are always degenerate. An earlier re- 
lated work based on Dirac theory considered both intrinsic and 
Rashba SOCs. The Vq dependence of G for the Rashba dom- 
inated case in that work agrees well with the MLG+R curve 
shown in Fig. 10, including the shoulder 



IV. EXPERIMENTAL ASPECTS 

A. Rashba spin splitting in graphene 

Whereas the Rashba spin splitting in MLG induced by an 
applied electric field is in general in the order of no more than 
100/ieV, which is beyond the present resolution of angle- 
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resolved photoelectron spectroscopy (ARPES), direct experi- 
mental observation of the Rashba spin splitting at K and K' in 
agreement with the first-principles calculations^"^ is so far not 
reported. An earlier experiment on epitaxial graphene layers 
on a Ni(l 1 1) surface reported a large Rashba interaction^*^ up 
to 225 me V but was soon questioned since the splitting might 
simply reveal a Zeeman type splitting due to the ferromagnetic 
nature of nickel.^' An intercalated Au monolayer between the 
graphene layer and the Ni(lll) substrate reduced the split- 
ting to about 13meV and was concluded as the Rashba effect 
on the Tt states supported by spin-resolved ARPES. How- 
ever, the low-energy band structure of MLGh-R at that time 
was not yet clear, and a simplified picture was adopted in the 
explanation of the measured spin splitting. In addition, trans- 
port properties of graphene based on metallic substrates can 
be difficult to isolate since a large bulk current will interfere 
as background. 

Throughout the above calculations we have mostly focused 
on a rather weak Rashba hopping parameter Ir = 30/ieV, 
yielding a splitting at the K and K' points 3fff = 90 /ieV, which 
is a realistic and rather conservative estimate for the gate- 
voltage-induced Rashba SOC strength. A recent proposal of 
impurity-induced SOC in graphene, however, indicated that 
the coupling strength can be strongly enhanced by putting 
heavy adatoms^^ as well as by hydrogenation.^^'^^ 

B. Klein tunneling in MLG 

Indirect and direct experimental evidences of Klein tunnel- 
ing in MLG have been reported recently.^^'^^. For detailed 
reviews, we refer to Refs. 2, 3, 20, 21, and 59. A very re- 
cent experiment on transport through a pnp junction in MLG 
used an embedded local gate, which yields high quality bal- 
listic transport and perfectly independent control of the local 
carrier density, as well as the feature of Klein tunneling. 

Recall the = curve of conductance for MLG shown in 
Fig. 10. Overall, the conductance for n — > n transmission with 
Vb < is always higher than that for n ^ p transmission with 
Vo > Ef. Even though Klein tunneling leads to perfect trans- 
mission at normal incidence in the latter case, the decay of 
T with incident angle eventually yields a lower conductance 
after integration. This feature has been agreed in recent ex- 
periments for pn and pnp junctions in MLG.'*"'^^'^**'^""^'* The 
difference of the conductance, or equivalently the resistance, 
between the nn and np (or between pp and pn) in experiments 
is even more obvious possibly due to the smooth interface that 
leads to an exponentially decaying form of T,^^ as we have re- 
viewed and discussed in Sec. Ill A 2. In fact, for MLG we have 
numerically checked G for pn junctions with a smooth inter- 
face, which indeed can enhance the difference of G between 
the nn and np regimes. 

Another interesting feature so far experimentally reported 
only in Refs. 58 and 60 is the Fabry-Perot oscillation of the 
conductance for pnp junctions due to the interference between 
the two interfaces of the central barrier. This feature requires 



the system to be ballistic and can be naturally revealed by our 
tight-binding transport calculation, which we will elaborate 
elsewhere in the future. 



V. CONCLUSION AND OUTLOOK 

In conclusion, we have employed tight-binding calcula- 
tions to show that transport properties of MLGh-R behave as 
BLG due to their identical form of the low-energy disper- 
sion, choosing the chiral tunneling in pn and pnp junctions 
as a concrete example. Within single-band transmission, nor- 
mal incidence transmission through a pn junction in BLG 
with \Ef\ < 7j is forbidden as a consequence of pseudospin 
conservation,^^ while in MLGh-R with \Ef\ < 3tR this forbid- 
den transmission also occurs but as a consequence of real-spin 
conservation. In mapping the angle- and spin-resolved trans- 
mission for the MLGh-R case, a combined effect of forbidden 
normal transmission and intrinsic spin-Hall deflection is re- 
vealed [Fig. 8(b)]. Compared to the potential-free spin-Hall 
deflection case as shown in Fig. 8(a), where = Ti = 1 at 
(/) = 0, the effect of the pn junction seems to force the up and 
down spins to separate since = = at = 0. The fea- 
ture revealed in Fig. 8(b) may therefore suggest a new type of 
intrinsic spin-Hall mechanism in MLG. 

Within multiband transmission, however, the Rashba SOC 
in MLG no longer plays an important role when \Ef \ ^ 3ffi 
(Fig. 7). Likewise, the interlayer hopping Yi in BLG becomes 
unimportant when \Ef | S> 7i . Transport in both MLG+R with 
I^fI > 3tR and BLG with l^f | > recovers to fliat in MLG, 
despite the usually very different energy scales of 3?^ and 7j . 
In view of the distinct transmission patterns in MLG+R with 
I^fI < 3tR [Fig. 8(b)] and \Ef\ > 3tR (Fig. 7), as an inter- 
esting conjecture for the BLG case one expects very different 
scattering regimes for \Ef \ < Yi and \Ef | S> 7i . The former is 
well discussed in the literature and exhibits strong scattering 
[Fig. 3(d)] while the latter is less discussed and the scattering 
is expected to be strongly suppressed. 

MLG and BLG are known to behave quite differently in 
general, in the sense of single-band transmission. Whereas 
turning MLG directly into BLG is in principle not possible, 
steering MLG to MLG+R can be achieved simply by gating, 
and therefore the effect of Rashba SOC provides a possibility 
to continuously change the MLG-like transport properties to 
BLG-like. We expect further transport properties to behave 
similarly in BLG and in MLG+R, such as the quantum Hall 
effect,^^ as was also noted by Rashba. 
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